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The unsteady flows caused by the interaction between the impeller and the volute in a high-speed micro centri- 


fugal pump are numerically studied. The internal flows of both with and without cavitations are analyzed using 


the CFX. The characteristics of unsteady pressure on the blade surfaces and the symmetric plane of the volute are 


presented and compared. The results show that the amplitudes of pressure fluctuations of critical cavitation on the 


blade pressure surface (PS) are bigger as compared with those at the non-cavitation condition, but on the suction 


surface (SS), the situation is on the contrary. When cavitation occurs, reduction of load in the impeller is a result. 


In the present study, such reduction of load is observed mainly on the first half of the blades. Pressure fluctuations 


at five monitoring points, denoted by WK1 to WKS in the volute, are also analyzed. No matter at the critical ca- 


vitation or at the non-cavitation conditions, the monitored pressure fluctuations are at the same frequencies, which 


equal to the blade passing frequency (BPF) and its multiples. However, the amplitudes of the fluctuations at criti- 


cal cavitation condition are considerably stronger, as compared with those for without cavitation. 
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Introduction 


Cavitation is a complex phenomenon of liquid vapo- 
rizing and forms bubbles in local low pressure regions. It 
is a common occurrence in hydraulic fluid machinery [1]. 
When cavitation occurs, the overall performance such as 
the efficiency of the pump may decline. Cavitations may 
also lead to some disadvantages such as vibration, noises 
and failure of impeller. In order to solve the extremely 
difficult problem relevant to the cavitation, the cavitating 
flow characteristics in a pump should be carefully studied. 

When a micro-centrifugal pump operates at high speed, 
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cavitation is very easy to occur in the flow passages of 
impeller. However, this kind of the internal flows in high 
speed micro-centrifugal pumps, especially the unsteady 
cavitation phenomena, are still very insufficiently studied 
as compared with those of traditional large scale centri- 
fugal pumps. 

There are many studies in which the cavitating flow 
fields in large scale centrifugal pumps are investigated 
[3-5]. For instances, Medvitz et al. [3] adopted a cavita- 
tion model that was proposed by Kunz et al.[2], to predict 
the internal cavitation flows in a centrifugal pump. The 
curve of head-drop was calculated at several flow rates. 
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Nomenclature 
b> blade width at exit[mm] T period of the impeller rotation[s] 
BPF blade passing frequency[Hz] Uy velocity at the impeller outlet[m/s] 
G, static pressure coefficient number of blades 
Dy diameenot hie impeller cule) dimensionless distance from the hub to the 
shroud 
Fond cavitation condensation coefficient Y, volume fraction of vapor 
Fiap cavitation vaporization coefficient N volume fraction of liquid 
H head P, density of vapor 
n impeller rotating speed[r/min] P, density of liquid[kg/m*] 
Pin static pressure of impeller inlet[Pa] n efficiency[%] 
Pins the static pressure root-mean-square Oo cavitation number 
Q flow rate[m*/h] w head coefficient 
Ou design flow rate[m*/h] © surface tension coefficient[N/m] 
s non-dimensional arc length at the midspan of the Ory, cavitation number at which the pump head 


blade, from the leading to the trailing edges 


Fu et al. [4-5] applied a cavitation model, which had been 
implemented in commercial software, to numerically 
study cavitation flow phenomena in a centrifugal pump 
under low flow rates. They observed that cavitation oc- 
curred over a wide range of low flow rates, which showed 
a typical head-drop curve. Axially asymmetric cavities 
distributed in the impeller passages were also observed. 

It is known that the rotational speed and the size of 
micro-centrifugal pumps are much different from those 
of conventional pumps. So the cavitation curves calcu- 
lated by the similarity law for ordinary large scale pumps 
are not suitable for a micro pump [1]. Compared with 
large scale low-speed pumps, a high speed micro-centri- 
fugal pump has a higher specific speed, which affects the 
cavitation significantly. In the present study, the internal 
unsteady cavitating flow characteristics of a high speed 
micro centrifugal pump are numerically studied by using 
a commercial package, the CFX. For the closure of the 
governing equations, the RNG k-e model for turbulence 
and the Rayleigh-Plesset model for cavitation are em- 
ployed. This study aims to analyze the influence of cavi- 
tation on unsteady pressure fluctuations in the high speed 
micro centrifugal pump. The flow details on blade sur- 
faces and in the symmetric plane of the volute will be 
paid attention. 


Model studied and numerical scheme 


Basic geometric parameters of the pump 


The computational model of the pump includes a 
shrouded impeller which has three blades, a volute, and 
the inlet and outlet pipes. The main relevant geometric 
and design parameters of the high speed micro centrifug- 
al pump are summarized in Tablel. 


is degraded by 3% because of cavitation 


Mesh and check of grid-independency 


Unstructured tetrahedral cells were used for modelling 
the twisted blades and the complex internal flow passag- 
es of the micro-centrifugal pump, because of their good 
extensive adaptability to complicated geometries and the 
easiness for grid generation. 


Table 1 Design parameters of the centrifugal pump 


Flow rate: Qa (m*/h) 13.2 
Head: H (m) 30 
Rotational speed: n (r/min) 10000 
Specific speed: ns 169 
Impeller diameter: Dz (mm) 56 
Blade outlet width: b2 (mm) 4 
Blade outlet angle: 52 (deg) 39 
Blade Number: Z 3 


(a) impeller 


(b) tongue of volute 
Fig. 1 Grids of the impeller and tongue 


The grids near the blade surfaces and in nearby of the 
volute tongue are refined to improve the resolution of the 
results. Figurel shows the refined grids of the impeller 
and the tongue. 
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Four meshes listed in Table 2 have been used to check 
the grid-independence of the calculation. The check is 
based on calculating the steady flow with the cavitation 
number equal to 0.115. General performance such as the 
head and the efficiency of the pump calculated by the 
four meshes are included in Table 2, which shows very 
slight differences between the results of meshes III and 
IV. Figure 2 is the comparison of the vapour fraction 
calculated by the four meshes, where S is the non-dimen- 
sional arc length of the midspan, from the leading to the 
trailing edges of the blade; Z, is defined as the dimen- 
sionless distance from the hub to the shroud, and its val- 
ue is in the range from 0 to 1. Moreover, as clearly seen 
in Figure 2, the calculated vapour volume fraction is re- 
mained constant when the grids are refined from mesh HI 
to mesh IV. Therefore, the mesh III which has 1682417 
elements is employed for the present simulations. 


Table 2 Meshes for grid-independence check 


Meshes 
zones 
I Ul lll VI 
Inlet 196918 196918 196918 196918 
Impeller 196400 369994 474467 603025 
Volute 778216 778216 778216 778216 
Outlet 232816 232816 232816 232816 
Total 1404350 1577944 1682417 1810975 
H/m 30.73 30.64 30.36 30.39 
n/% 76.15 76.82 71.32 77.27 
í 
——— =: Mesh | 
O08 eR ste Mesh II 
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Fig. 2 Calculated vapor fraction along the line Z,=0.5 in SS of 
the blade 


Cavitation model 


The cavitation flows of the fluid are assumed to be a 
homogeneous two-phase mixture in the cavitation model. 
The liquid-vapor mass transfer due to cavitation is solved 
by the transport equation-based cavitation model. 

AY, Py) 
ôt 

Where % is the volume fraction of vapour, p, is the 


+ V(7,p,u) =m, — m (1) 
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vapour density, m, and m, are the mass transfer rates 


per unit volume for vaporization and condensation, re- 
spectively. The values of the m, and m, are deduced 
from the Rayleigh-Plesset n 


2 = 
Ry u e P yy 
d 2 dt 7 Pı 
where R, is the bubble radius, P is pressure while 


Psat is the saturation pressure for the vapour, p; is the liq- 
uid density and © is the surface tension coefficient. By 
ignoring the effect of the second order derivative and the 
surface tension, the Rayleigh — is simplified as: 
dR, 2 Pu- 
Pa (Eiu iy (3) 
dt 3 P/i 
Denote the number of bubbles per unit volume by Nz, 
then the rate of bubble mass transfer pe unit fona is: 
dm d Vg 1/2 
2 =N - Npp, 4R (E a P (4 
a a a ae aG a fsa") (4) 
The values of Ng for vaporization and condensation, 
respectively, are calculated as follows [6,7], 


3a 
Np =e ~g) (5) 
a 4nR; 
3a 
Np = 4 6 
a 4r R? (9) 


The mass transfer rates per unit volume for vaporiza- 
tion and condensation, respectively, are therefore calcu- 
lated by 


Epp DP Ela EN Pepu (N) 
P Rg PI 
3 2 P-P 
m, = Peona da A sar e P> Eat (8) 
Rpg PI 


In the Rayleigh-Plesset cavitation model, me = 5 x 
10%. Fap and Fong are empirical constants, Fap = 50, and 
Fond = 0.01, as recommended by Zwart et al.[6] and 
Bakir et al. [7]. 


Computational scheme 


The unsteady cavitation flow simulations are per- 
formed by using the commercial software CFX. The 
RNG k-e model for turbulence and the Rayleigh-Plesset 
model for cavitation are selected for calculation. The 
RNG k-e turbulence model is widely used to account for 
turbulence in cavitating flows [8-11]. The vapour and 
liquid volume fractions are zero and one, respectively, at 
the inflow boundary the total pressure is also specified. 
Target mass flow is applied at the outlet of the computa- 
tional domain. In addition, no-slip and adiabatic condi- 
tions have been imposed on the solid walls. The transient 
rotor/ stator mode is applied at the interfaces between the 
impeller and volute. 180 time steps are set in the period 
of a rotation T; this results in a time step for unsteady 
calculation, ts=T/180=3.333*10° (s), which is equivalent 


LI Yeqiang etal. Numerical Study of Unsteady Flows with Cavitation in a High-Speed Micro Centrifugal Pump 21 


to the time interval for the impeller to rotates 2°. In order 
to ensure the convergence of calculation, the temporal 
history at a monitoring point which locates in the inter- 
section between the midspan rotary surface (Z,=0.5) and 
the impeller/ volute interface is recorded and shown in 
Figure 3. The static pressure coefficient fluctuation at the 
monitored point has shown periodical unsteady changes 
after ten cycles (1800 iteration steps), which indicates 
that the flow is fully developed and the calculation is 
converged. The static pressure coefficient is defined as: 


P-P 
Cp =— 5 (9) 
PiU2 / 2 
0.75 
0.7 
0.65 
S 
0.6 
0.55 
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Fig. 3 Unsteady pressure signal at monitoring point 
Results and Discussion 


Cavitation performance 


The cavitation performance of the micro-centrifugal 
pump is experimentally tested according to the Chinese 
GB3216-82 Standard. The test rig is a closed circuit with 
a container in the suction side of the pump, while a va- 
cuum pump is installed on the container so as to adjust 
the suction pressure of the tested pump. The flow-rate is 
measured by an electromagnetic flow meter (BFG- 
SDNSO type, range from 0 to 20 m*h ', an accuracy over 
the calibration range of 0.18% of reading). The pressure 
head is measured using a high resolution Bourdon tube. 
Input shaft-work is measured by calculating the 
work-output of the driving motor, whose rotation speed is 
monitored and controlled at about 10000 r/min. 

Characteristics of cavitation in the high-speed micro 
centrifugal pump are simulated at several flow-rates. 
Figure 4(a) shows coefficient degradation characteristics 
of the head when cavitations occur, at 80%, 100% and 
120% of the design flow rate and under different cavita- 
tion numbers. The cavitation number and head coefficient 
are defined as: 


pau E 
Pı 
=s (11) 
uy 


It can be observed that the head coefficient has a sud- 
den decrease when the cavitation number is decreased to 
the critical value, as expected under cavitating conditions. 
The critical cavitation number 03, which denotes the 
cavitation number at which the pump head is degraded 
by 3% because of cavitation, is obtained by simulation 
and experiment, and are shown in Figure 4(b). The nu- 
merical values of oy, are generally in good agreements 
with the experimental data. The results demonstrated that 
the numerical model and method can predict the cavita- 
tion flows precisely. 


Analysis of unsteady cavitating flow 


Figure 5 shows the distribution of volume fraction of 
vapour at the midspan of the impeller passage, for design 
flow-rate but with different cavitation number values. 


0.65 


06 PO 
=: 
0.55 
E —yw— 120d 


0.1 
—A— o CFD 
0.08 —E— o EXP 
0.06 
x 
6 
0.04 
0.02 


0.8 0.9 


1 1.1 2 
Q/Qd 
(b) 0-63», 


Fig. 4 Performances of the micro centrifugal pump 
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Fig.5 Vapour fraction in the Z,=0.5 of the impeller 
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A small and thin vapor bubble is initially observed to 
form near the leading edge of the suction side of the 
blades, as can be seen in Figure 5(a). As the cavitation 
number drops to 0.054, the vapor bubbles grow in size 
along the suction side of the blades, and stretch to the 
middle of the blades in the flow passages. When the cavita- 
tion number is further decreased to 0.045, which is called 
critical cavitation condition under design flow rate, the 
vapor bubbles almost occupy a half of the flow passages 
between blades, while the flow in the impeller is seriously 
affected. As a result, the head of the pump is degraded. 
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Figure 6 shows the instantaneous pressure distribution 
in the critical cavitation condition, in both the PS and SS 
sides of the blade and at three spanwise locations. Seeing 
from Figure 6, local minima appear on pressure surface 
of the blade near the leading edge. The pressure increases 
along the blade firstly, than start to decrease around 
S=0.8. Meanwhile, with the time changing, the pressure 
on blade surface will increase gradually; this is due to the 
blade moves close to the volute tongue. The pressure on 
blade suction surface at different Z, locations is remained 
low values of almost 0. 


PS,t=27/3 4", 
SS,t=27/3 ," 
Po 


Fig. 6 Instantaneous pressure profiles along blade surface 


In order to analyse the effects of critical cavitation on 
the load of blade, Figure 7 shows the distribution of the 
blade surface loads at the midspan, at both non-cavitation 
and critical cavitation conditions. The blade surface load 
is defined as the pressure difference between the pressure 
and the suction sides at the same radius. Figure 7 clearly 
shows that, at both of the two conditions, the blade loads 
display a trend of increasing firstly and then decreasing, 
from the blade leading edge to the trailing edge. Howev- 
er, in the first half of the blade, the blade load C, at the 
critical cavitation condition is considerably lower than 
that of the non-cavitation case. So when the cavitation 
occurs and develops, the work-doing capacity of the 
blades in the first half of their length is cut down. 

Figure 8 shows the distribution of the root-mean- 
square of the unsteady pressure, the P,,,, on the blade 
surface. P,„ms 1s defined as follows: 


N E 
DE - P? jn 
pad 


rms pus /2 

where N is the number of time steps in a rotation of 

the impeller, N=180 in the current calculation as men- 
tioned above. 

It can be seen in Figures 8 (a) and 8(c) that P,ms grad- 

ually decreases from the trailing edge to the leading edge 

on blade surfaces, while along the blade height direction 


(12) 
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Fig. 7 Blade load along midspan arc-length 


P,ms Changes very slightly. The P,ms denotes the circum- 
ferential non-uniform pressure fluctuations, caused by 
the interaction between the impeller and the volute. The 
P,ms Value of critical cavitation is larger than that at the 
non-cavitation condition, as shown in Figures 8(a) and 
8(b). However, on the suction surface, the situation is on 
the contrary. As shown in Figure 8(d), the P,ms is almost 
zero, except at the trailing edge of the blade for the criti- 
cal cavitation case. In such situation, due to the genera- 
tion and development of vapour on the suction surface of 
blade, the pressure fluctuations become weak within the 
vapour region, the value of P,ms is therefore small. 
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0.206 114182226 3 


(a) PS at non-cavitation 


0.206 1 14182226 3 


(b) PS at critical cavitation 


0.20.6 1 1.41.8 2.2 


(c) SS at non-cavitation 


0.20.6 1 1.41.8 2.2 


(d) SS at critical cavitation 


Fig.8 Pns on PS and SS of the blade 


In order to analyse the non-uniform circumferential 
fluctuations, 5 monitoring points are set up in the sym- 
metric plane of the volute; positions of the monitoring 
points are shown in Figure 9. The fast Fourier transform 
(FFT) results for the frequency characteristics of pressure 
fluctuations at design flow rate for both the non-cavita- 
tion and critical cavitation flows are given in Figure 10. 
The dominant frequency of the pressure fluctuation is the 
blade passing frequency (BPF) at all monitored points for 
both the non-cavitation and the critical cavitation flows. 
The fluctuation amplitudes at the critical cavitation con- 
dition are considerably stronger than those of without 
cavitation case. Comparing the pressure modes in the 
frequency domain, a peak value is firstly found at BPF at 
monitoring point WK1 which is near the tongue of the 
volute. The maximum amplitude of the pressure pulsa- 
tion is found at monitoring point WK2, which is about 50 
degree in downstream of the tongue. 

In order to find the reason why pressure fluctuation at 
WK2 is stronger than that at WK1, secondary flow 
streamlines distribution at WK1 and WK2 sections in the 
critical cavitation case are compared in Figure 11. The 
results show that the secondary flows at WK1 are weaker 
than those at WK2, where large-scale vortices are ob- 
served all the time. Reasonably, the stronger instantane- 
ous secondary motions in cross-section WK2 contribute 
to the intense pressure fluctuations. 


Fig.9 Monitoring points 
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(b) critical cavitation 


Fig. 10 FFT results of static pressure at monitoring points 
Conclusions 


In this paper, numerical study on the unsteady cavita- 
tion flow in a high speed micro centrifugal pump has 
been carried out. The main conclusions are summarized 
as follows: 

The amplitudes of pressure fluctuation at critical ca- 
vitation condition on the blade pressure surface are 
stronger than that in the non-cavitation case; but on the 
suction surface of the blade, the situation is on the con- 
trary. When the cavitation occurs, reduction of load is 
observed mainly on the first half of the blades in the 
flow passages. 
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(b) WK2 section 


Fig.11 Secondary flow in WK1 and WK2 sections 


The unsteady interaction between impeller and volute 
results in large scale pressure fluctuations. When cavita- 
tion occurs, such interactions become stronger in ampli- 
tudes. Peak pressure fluctuations due to unsteady interac- 
tion are firstly observed near the volute tongue, while the 
strongest fluctuations occur at about 50 degrees in down- 
stream of the tongue where the secondary flow is also 
strong and unsteady. 


Acknowledgments 
The study in this paper is funded by the National Nat- 


ural Science Foundation of China under Grants 
51176048 and 51576067. 


J. Therm. Sci., Vol.26, No.1, 2017 


References 


[1] X. Guan. Modern Pumps Theory and Design. China As- 
tronautic Publishing House, Beijing, 2011: pp.74-80. 

[2] R. F. Kunz, D. A. Boger, D. R. Stinebring, et al. A pre- 
conditioned Navier-Stokes method for two-phase flows 
with application to cavitation prediction. Comput & flu- 
ids 2000; 29(8): 849-875. 

[3] R. B. Medvitz, R. F. Kunz, D. A. Boger, et al. Perfor- 
mance analysis of cavitating flow in centrifugal pumps 
using multiphase CFD. Journal of Fluids Engineering, 
2002, 124(2): 377-383. 

[4] Y. Fu, J. Yuan, S. Yuan, et al. Numerical and Experi- 
mental Analysis of Flow Phenomena in a Centrifugal 
Pump Operating Under Low Flow Rates. Journal of Flu- 
ids Engineering, 2015, 137(1): 011102. 

[5] J. Yuan, Y. Fu., S. Yuan, A Study of Cavitation Flow 
in a Centrifugal Pump at Part Load Conditions Based on 
Numerical Analysis. ASME 2012 Fluids Engineering 
Division Summer Meeting collocated with the ASME 
2012 Heat Transfer Summer Conference and the ASME 
2012 10th International Conference on Nanochannels, 
Microchannels, and Minichannels. American Society of 
Mechanical Engineers, Rio Grande, Puerto Rico, 2012: 
193- 202. 

[6] Zwart P. J.,Gerber A. G. Belami T. “2004 A two phase 
flow model for predicting cavitation dynamics,” In: Pro- 
ceeding of the 5th International Conference on Maulit- 
phase Flow (ICMF 2004). 30 May~ 3 June 2004, paper 
no. 152. Yokohama, Japan: Japanese Society for Multi- 
phase Flow. 

[7] F. Bakir, R. Rey, A. G. Gerber, et al. Numerical and ex- 
perimental investigations of the cavitating behavior of an 
inducer. International Journal of Rotating Machinery, 
2004, 10(1): 15-25. 

[8] V. Yakhot, S. A .Orszag, S. Thangam, et al. Develop- 
ment of turbulence models for shear flows by a double 
expansion technique. Physics of Fluids A: Fluid Dynam- 
ics, 1992, 4(7): 1510-1520. 

[9] L. Zhou, Z. Wang. Numerical Simulation of Cavitation 
Around a Hydrofoil and Evaluation of a RNG x-e Model. 
Journal of Fluids Engineering, 2008, 130(1): 011302. 

[10] L. Tan, S. Cao, Y. Wang, et al. Direct and inverse itera- 
tive design method for centrifugal pump impellers. Pro- 
ceedings of the Institution of Mechanical Engineers, Part 
A: Journal of Power and Energy, 2012, 226(6): 764-775. 

[11] N. Dittakavi, A. Chunekar, S. Frankel. Large Eddy Si- 
mulation of Turbulent-Cavitation Interactions in a Ven- 
turi Nozzle. Journal of Fluids Engineering, 2010, 132(12): 
121301. 


